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Abstract 

Motivated by the fluid/gravity correspondence, we consider the Penrose inequality in the frame- 
work of fluid dynamics. In general relativity, the Penrose inequality relates the mass and the entropy 
associated with a gravitational background. If the inequality is violated by some Cauchy data, it 
suggests a creation of a naked singularity, thus providing means to consider the cosmic censorship 
hypothesis. The analogous inequality in the context of fluid dynamics can provide a valuable tool 
in the study of finite-time blowups in hydrodynamics. We derive the inequality for relativistic and 
nonrelativistic fluid flows in general dimension. We show that the inequality is always satisfied at 
the ideal fluid order. At the leading viscous order, the inequality may be violated by relativistic 
fluid flows, while it is always satisfied by nonrelativistic incompressible flows. The inequality may 
be violated at the next to leading viscous order by both relativistic and nonrelativistic flows. 

PACS numbers: 04.70.-s, 47.10. ad, 11.25.Tq 
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I. INTRODUCTION 



Compressible hydrodynamics allows singularities in the form of shock waves, that is a 
jump in the velocity vector field. Incompressible hydrodynamics does not exhibit shock 
waves, and by singularities one means blowups of the derivatives of the velocity vector field. 
It is often expected that the viscosity terms act as regulators that prevent the blowups. 

The basic question concerning singularities in the hydrodynamic description, is whether 
starting with appropriate initial conditions, where the velocity vector field and its derivatives 
are bounded, can the system evolve such that it will exhibit within a finite time a blowup 
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of the derivatives of the vector field. Mathematically, this important question has been 
presented as one of the millennium problems posed by the Clay Mathematics Institute 
Physically, such singularities if present, indicate a breakdown of the effective hydrodynamic 
description at long distances and imply that some new degrees of freedom are required. 

The issue of hydrodynamic singularities has an analogue in gravity. Given an appropriate 
Cauchy data, will the evolving space-time geometry exhibit a naked singularity, i.e. a blowup 
of curvature invariants and the energy density of matter fields at a point not covered by a 
horizon. The cosmic censorship conjecture in general relativity (for a review see [2J), states 
in its weak form that any space-time singularity is generically enclosed by an event horizon, 
that is a region that no light rays can escape from. Mathematically, naked singularities can 
exist, but physical considerations of causality seem to require that it should not be visible 
to distant observers. In its strong form, the conjecture states that even inside the black 
hole a falling observer will never "see" the singularity, i.e. time-like singularities never occur 
generically. 

n n 

The Penrose inequality [3] (for a review see e.g. [4j) is a conjecture relating the mass 
and the horizon area of any Cauchy initial data that if violated leads to a space-time naked 
singularity. Let us briefly review Penrose's argument in asymptotically flat four-dimensional 
space-time. Consider a Cauchy data (.M 3 , g, K), where A4 3 is a 3-manifold with a complete 
Riemannian metric g and an extrinsic curvature K, satisfying the Gauss-Codazzi equations. 
In particular, let an asymptotically flat initial Cauchy data with a mass Mo and event 
horizon area Aho evolve with time. The solution to Einstein equations with this initial data 
is expected to settle down at long time to a Kerr black hole solution with mass M and 
horizon area A H , where the inequality (we use c = Gjy = 1) 

M > ^/Ah/16% (1) 

holds. The inequality is saturated by the Schwarzshcild black hole. 

By the Hawking area theorem 5|, the event horizon area does not decrease with time 
Ah > Aho- The mass, on the other hand, cannot increase and may only decrease due to 
radiation loss M < Mq. Thus, we get the chain of inequalities 



M > M > ^A H /I6n > y/A H o/lGir , (2) 



and the initial Cauchy data should also satisfy the Penrose inequality M > a/ Ahq/IQk- 



The argument relies on the Hawking area theorem and the relaxation at late times to a 
Kerr solution, both assume the weak censorship hypothesis. Therefore, finding Cauchy data 
that violates the Penrose inequality implies finding a solution to Einstein equations in which 
a naked singularity is created. Proving that the inequality is satisfied does not prove the 
cosmic censorship conjecture but provides an important support for it. Note, that there are 
other less likely ways to interpret a violation of the inequality, e.g. that the system never 
reaches equilibrium. 

In fact, we can formulate an analog of the Penrose inequality for any dynamical system. 
Consider a system, where the total energy does not increase. We replace the area increase 
theorem by the second law of thermodynamics, that is the total entropy never decreases. 
The system is assumed to relax at late time to a thermodynamical equilibrium state. The 
inequality relation between the mass and the entropy is simply an equilibrium relation, while 
by repeating Penrose's argument, the inequality should hold away from equilibrium. 

The fluid/gravity correspondence relates in a derivative expansion two non-linear systems, 
the fluid dynamics described by the Navier-Stokes equations and gravity described by general 
relativity. In one particular (conformal) setup of the correspondence (J, the Navier-Stokes 
equations in d space-time dimensions are identified with the constraint Einstein equations 
of a deformed black brane gravitational background in one higher dimension. Alternatively, 
they are identified with the Gauss-Codazzi equations governing the evolution of the black 
brane event horizon j3"9|. Motivated by the fluid/gravity correspondence, we consider the 
Penrose inequality in the framework of fluid dynamics, which can provide a valuable tool in 
the study of finite time blowups in hydrodynamics. We derive the inequality for relativistic 
and nonrelativistic fluid flows in general dimension and study its implications. 

The paper is organized as follows. In section II we briefly review the fluid/gravity corre- 
spondence in the conformal case. We then derive the inequality for relativistic ideal hydro- 
dynamics, where we consider both the conformal and non-conformal cases. We show that 
the inequality is always satisfied at the ideal fluid order, and present a class of solutions 
that develop a shock wave singularity, while still satisfying the inequality. This is curious, 
though of course there is no inconsistency, since a non-violation of the inequality is only a 
sufficient condition for not having a finite-time blowup and not a necessary one. In section 
III we derive the inequality at the first and second viscous order for relativistic fluids. The 
inequality is nontrivial and maybe be violated. In section IV we show that the inequality is 
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always satisfied by nonrelativistic incompressible fluids. However, this is not the case when 
considering subleading corrections to the incompressible Navier-Stokes equations. Section 
V is devoted to a discussion. 



II. THE PENROSE INEQUALITY IN HYDRODYNAMICS 

In this section we will formulate the Penrose inequality for fluid dynamics. We will use 
the fluid/gravity correspondence to motivate the inequality. However, the actual derivation 
does not require the use of the correspondence. 



A. The fluid/gravity correspondence: conformal case 

Consider the vacuum (d + l)-dimensional Einstein equations with a negative cosmological 
constant 

RMN + dg MN = 0, R=-d(d+l) (M, N — 0, .., d) . (3) 

We will use the notation X M = ji = 0, d — 1, x d = r). One solution of these equations 
is the boosted black brane with AdSd+i asymptotics 

dsl = —lu^dx^dr + d2 u tl u u dx^dx u + r 2 rj^ v dx^dx u , (4) 

with the d- dimensional velocity vector u M = (7, JV{), u^u^ 1 = — 1, b related to the Hawking 
temperature by 

b = -4= , (5) 

and r]^ = diag(—l, 1.., 1) is the d- dimensional Minkowski metric. 

The metric (j4]) is an exact solution of the equations (J3j) when b and are constant. 
This is the thermal equilibrium state. Next, we get out of equilibrium by deforming the 
solution and allowing With these, the metric (BJ is no longer a solution of the 

equations (j3J). 

One can, however, correct the metric 01]) by adding derivatives of We will 

have in general 

ds 2 = dsl + dsl + ds\... (6) 
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where the lower index of ds 2 denotes the order of the derivative terms that is contains and 
dsl is the metric (j3J. Requiring the deformed metric (JBJ) to be a solution of (J3J) up to 
the neglected order in derivatives, is equivalent the order by order derivative expansion of 



conformal hydrodynamics jfj]. More specifically, the set of d+1 constraint Einstein equations 
can be written as the conservation law of a conformal hydrodynamic stress energy tensor in 
flat d- dimensional space with specific transport coefficients 

dfjT^ = . (7) 

Thus, the gravity solution provides a dual description in terms of fluid dynamics. Equiva- 
lently, one arrives at the fluid/gravity correspondence starting from the AdS/CFT correspon- 
dence, and considering the hydrodynamics of the CFT and its gravitational dual description. 
Alternatively, the hydrodynamic equations are identified with the Gauss-Codazzi equations 



governing the evolution of the black brane nul 
related by an RG flow in the radial direction 



event horizon 



io|. 



7H9| . The various schemes are 



The hydrodynamic stress energy tensor TL, can also be calculated at the <i-dimensional 



boundary 11] . We write the (d + l)-dimensional metric in an ADM-like decomposition 

ds 2 = N 2 dr 2 + ^(dx* + N^dr)(dx u + N v dr) , (8) 

where 7 VA1 is the boundary (r = const.) metric, N and are the lapse and shift functions. 
The extrinsic curvature tensor of the boundary DM is defined by 

= -—{d rl ^ - V^N V - V„Ay , (9) 

and the boundary stress energy tensor reads 

r d-2 

Tfj, u = lim r -+ 00 [-^-2(K flv - Kj^ - (d - ■ (10) 

B. The Penrose inequality in conformal hydrodynamics 

The hydrodynamic stress energy tensor of an ideal relativistic neutral fluid is 

Tfu, = (e + p)UfjU u + prj^u , (11) 

where e is the energy density and p is the pressure. Conformal invariance implies that the 
stress energy tensor is traceless, thus we have the equation of state e = (d — l)p. The 
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only dimensionfull parameter is the temperature T and we have that e = (d — l)p ~ T d . 
The boundary stress tensor f lTUj) corresponding to the stationary metric (j3J) is that of ideal 



conformal neutral fluid and reads 
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Tp, = 7^p(^ + • ( 12 ) 

Let us derive the analog of the Penrose inequality in this case. To this end we need 
to define the energy (mass), the horizon area and the equilibrium relation. The energy 
(E = M) is given by 



E = j T 00 d^x = J d ^ x , (13) 

while the horizon area is defined via the entropy S. Consider a conformal fluid at rest 
(Ufj, = (1,0,0,0)) satisfying the thermodynamic relation 

e + p = Ts, (14) 

where s is entropy density, which in the conformal case reads 

6 d 1 (15) 



Td-1 AW- 1 
and the horizon area is 

1 



A = AS = J ^d^x , (16) 

so we get that the equilibrium relation between energy and entropy 

E = d ^S)^ , (17) 

where we normalized the volume of space to one. The fluid analog of the Penrose inequality 
is therefore 

E>^(4S)^. (18) 

One way to verify the direction of the inequality is to consider the black brane with angular 
momentum. 

Consider as an example the case when the fluid is not at rest. The inequality is slightly 

generalized in this case by replacing the entropy density s by the time component of entropy 

^ fl4J ]. that is s° = 57. For any 7 > 1 we have 

g?7 2 —Id . . 

\ : > 7^ • (19) 



Thus, by Holder's inequality we get 

f ^ ~ l > I ^ > ( / ^-)^ • (20) 
J b d (d-l) 'J ¥ - V b d - l> { ' 

The inequality (|T8|) holds for any 7 > 1, and is an equality only when 7 = 1, that is for a 
fluid at rest. 

Note, that although motivated by the fluid/gravity correspondence the argument that 
leads to the inequality for fluids is rather general and is independent of having a gravitational 
description. It assumes that at late times the system relaxes to an equilibrium state, whose 
thermodynamic relation gives the inequality relation. One then assumes that mass (energy) 
is a non-increasing, that is no energy is pumped into the system, and that the entropy 
is non-decreasing by the second law of thermodynamics. Thus, starting with initial data 
(E , So) one has 

E > E > ^(45)^ > ^-1(AS )^ . (21) 

107T 107T 

The Penrose inequality is a condition to be imposed on any initial spacetime configuration 
(Cauchy data). When this condition fails, the evolution in time is expected to exhibit a 
naked singularity. Likewise, The inequality f fl8|) is a condition that we impose on any initial 
fluid data that satisfy the hydrodynamic equations. If the condition fails, we expect the fluid 
to exhibit a finite-time blowup of the velocity field gradients. Note, however, that a violation 
of the Penrose inequality is a sufficient condition for a creation of a naked singularity, but is 
not a necessary one. Likewise, a violation of the hydrodynamic inequality is only a sufficient 
condition for a singularity. 

It is quite important not to link directly singularities in the hydrodynamic description 
to naked singularities in gravity. The reason being that the fluid/gravity correspondence is 
established order by order in a derivatives expansion. An appearance of singularities signals 
the breakdown of a derivative expansion and for a direct link between the two descriptions 
one needs to go beyond the established relation. 

C. Nonconformal fluids 

We can easily generalize the inequality to non-conformal fluids with a general equation of 
state. As an example consider an equation of state of the form E = ap where a ^ d — 1. An 
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example of a dual gravitational background for such an equation of state was constructed 



in 



15] . A study of the universal statistics of such hydrodynamics in d = 2 has been carried 



out in 



161 ] . By dimensional analysis p = aT and 



e = oaT\ s = = (a + l)aT d - 1 . (22) 



T 

The relation between the energy and entropy reads 



O d , 

E = ~^^S^ , (23) 

(a + l) d -!a d -i 

and the inequality in this case is 

E > ^^S^ . (24) 

The inequality ffT8j) is a special case of ( ]2"4"|) with a = d — 1 and a = (^f-) d . 



D. Simple relativistic waves 

As stated above, a violation of the Penrose inequality is a sufficient condition for a creation 
of a naked singularity, but is not a necessary one. Likewise, a violation of the hydrodynamic 
inequality is only a sufficient condition for a singularity and not a necessary one. In this 
subsection we give an example of initial relativistic conformal ideal fluid flow configurations 
leading to a formation of a singularity, despite satisfying the Penrose inequality. 

The hydrodynamic equations for a relativistic conformal ideal fluid can be written in the 
form 

d^u" + {d- \)DlnT = 0, a a + P^d^lnT = , (25) 

where D = u^d^, a a = u fl d^u a . One obtains this form of the equations by the projection of 
the stress energy conservation equations as u u d^T tJ,u and P£d fl T tJ,u , with 

P(W = VfiU + U^Uy . (26) 

Consider the simple relativistic waves presented in [17] (here we generalize to d space-time 
dimensions) 



v + Cs n , fl + v 



u" = (r/(t,z),0...,0,>yv(t,z)), v(t,z) = f(z-—^t), b = — - ) . (27) 

1 T" VC S \ 1 V 
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where b is related to the temperature T as in (J5J), / is an arbitrary function specifying the 
initial wave profile, and c s is the speed of sound 

This configuration is a solution to equations f l25|) . and it disperses due its implicit de- 
pendence on v, which causes any initial regular single- valued wave to become ultimately 
multiple- valued, corresponding to the crossing of characteristics [171 ] . Therefore it develops 
a singularity for any regular initial conditions function f(z). The inequality (fT8|) is satisfied 
by the same argument as ( |20|) . and gives no indication of the singularity that will develop 
in the future. 



III. THE VISCOUS ORDER 

In this section we will derive the hydrodynamic inequality at the first and second viscous 
order for relativistic fluids. The inequality that we will construct is nontrivial and can 
potentially be violated. 



A. first order in derivatives 

The correction to the relativistic hydrodynamic ideal stress-energy tensor ffTTj) at the first 
viscous order reads 

r« = -2t?<v - (d a u a { V ^ + u M u„) , (29) 

where we work in the Landau frame u^T^J = 0. a^ v is the shear tensor 

11 1 
<V = -{dpUv + d v Ufj) + -{u^u, a d a u u + u v u a d a u^) - -j—^d a u a (r]^ u + u^u v ) . (30) 

7] > and £ > are the shear and bulk viscosities, respectively. In the conformal case, 
Tjf = implies that ( = 0. For the black brane gravitational background, the stress-energy 
tensor is 

Tilv = 16^ ( ^ + dUilUv) ~ 8^ a(lu ' (31) 
and the ratio of the shear viscosity to the entropy density is 2 = i. In the following we 
consider the conformal case. 
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The energy decreases due to shear effects and reads 



E= I [^n^-2 W oo) d d ~ l x, (32) 



16irb d 

where the first and second terms are the ideal and viscous contributions, respectively. The 
decrease in the energy is written explicitly 

- J 277(7 Q0 d d ~ l x = -2 J V (d tl + 7 (-7<9t + v%)j - -j^(-dtl + d^x . (33) 

The entropy, on the other hand, does not change at this order, but is no longer conserved 
~ ifO^v®^ ■ On the gravitational side, the horizon area of the deformed black brane is 



not corrected at the first viscous order 



6|- 



Thus, at the first viscous order the inequality reads 

" ' " 1 d d ~ x x - 2 / 77 (d a + 7 (- 7 $ + tM) 7 - £^±(-dt>y + d iV l ) ) d d ~\ 



ItotV J ' V d 

> — ( f ITT ^^)^ - (34) 

and may be violated. Shock waves regularized by the viscosity term of the type discussed 
in [3], may be possible candidates for violating the inequality (151]) . 

Note, that a fluid flow configuration that violates (151]) will also imply by the fluid/gravity 
correspondence that the deformed black brane background at first order in the derivatives 
will exhibit a naked singularity. This implies on the gravity side a need for the inclusion of 
the higher order derivative terms, possibly an infinite number of them. 



B. Bjorken flow 

As a simple example one can check that the Bjorken flow (for its description see e.g. []j]]) 
satisfies the inequality at the viscous order. The Bjorken flow is a one-dimensional expansion 
of plasma with boost invariance symmetry along the expansion axis. It is convenient to 
parametrize the flow by the proper time r and the rapidity y related to the lab frame 
coordinates by 

a; = rCosh[y], x 1 = rSinh[y] . (35) 

In equilibrium the energy density and the entropy density are e = eoT d and s = -^toT d ~ l . 
The inequality can be read from (1241) with a — d — 1 and a = The solution of the 
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hydro dynamic equations is obtained as an expansion in r~ a ,a = It is straightforward 
to check that the inequality is satisfied at the subleading order in this expansion, which 
captures the first order viscous corrections. 

C. Second order in derivatives 

Consider next the inequality at the second viscous order. It is straightforward to consider 
the most general stress-energy tensor at this order. However for simplicity we consider the 
particular form of the stress-energy tensor corresponding to the deformed black brane grav- 
itational background, that is a conformal hydrodynamics with a particular set of transport 
coefficients. The energy reads 

E = I T7TIrf(^ 2 - 1 - 2b(r 00 - 2br UJ (u x D x a 00 + 2w V A0 ) 

J 1O7T0 

+2b 2 (u x D x a 00 + a V A0 - ( 7 2 - 1))^^)) d d ~ l x . (36) 



r w is a new transport coefficient that appears at this order, in the vorticity tensor 

£«V = 2^ Uv ~ duU v) + 7;(uvU a d a u u - u v u a d a u^) , (37) 
and is the Weyl-covariant derivative that acts on the velocity vector as 

dx = (38) 

Then the entropy current at this order reads [f| 

S = J -^ ZT h + b 2 1 (A 1 a a ^ + A 2 u al3 u^ + A 3 R) 



+ b'^B.D^ + B 2 D x u 0X )] d d ~ l x , (39) 



where 



with 



R = 2{d- 1)8^ - (d - 2)(d - l)A^ , (40) 



d u a 

A^ = u a d a u^ - -j—^ u v ■ ( 41 ) 
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The coefficients A's and B's are defined in jfj] (Eq. (2.8) and Eq. (5.16)). 
The inequality at second order in derivatives is 

/ 16^ ^ ~ 1 ~ 2ba °° ~ 2bT ^ uXD ^oo + 2tu x a X0 ) 



+ 2b 2 (u x D x a 00 + a x a xo - ( 7 2 - 1))^^)) d^x > 

d-1 



16/T 



(/ ftirlT + b^iA^a^ + A 2 u a ^ + A 3 R) 
+ b 2 (B 1 D x a 0X + B 2 D x u 0X )] d^x)^ . (42) 



IV. NONRELATIVISTIC FLOWS 

In this section we will consider the hydrodynamic inequality for nonrelativistic flows. We 
will show that the inequality is always satisfied by nonrelativistic incompressible fluids, but 
not when considering subleading corrections to the incompressible Navier-Stokes equations. 



A. Navier-Stokes equations 

The incompressible Navier-Stokes equations can be obtained in the nonrelativistic limit of 
relativistic hydrodynamics, where the fluid flow velocity «<Cc. Specifically we can consider 
relativistic CFT hydrodynamics with the stress-energy tensor (l3Tj) . The speed of sound in 
relativistic CFT hydrodynamics v s = and thus, we are considering in the nonrelativistic 
limit also the low mach limit of velocities much smaller than the speed of sound. This is the 
reason for obtaining incompressible flows. 

The conservation law equations of the stress-energy tensor (IHTj) . projected on u v and P° 
read 



OX + (d- l)DlnT 



1 



2ttT 



0>c 



a a + P^lnT = _LpM(0 oO « _ ( d _ i )a « aa) . ( 43 ) 

In the nonrelativistic we expand w M = (1 — v 2 /2 + v *), T = T (l + P + ...) where we scale 
v l ~ e, di ~ e,d t ~ e 2 ,P ~ e 2 , where e ~ 1/c ~ l/v s . The first equation in ( 14"3]) gives the 
incompressibility condition diV % = 0, while the second equation gives [20I, |2l| 

d t v l + tfdjV* = ~d l P + i/Au* , (44) 
13 



where v = 



B. First order inequality 

The NS equations (jUJ) are of order e 3 . The inequality at the first viscous order is of order 
e 2 . Let us derive it. We use b = b (l — P), b = j4p- = dv. The energy up to order e 2 reads 

^^/O + ^+^K 1 *' (45) 

where we put c = 1. The entropy up to order e 2 reads 

AS = -L J (i + V l + (d - l)p) d d - l x . (46) 
It is easy to see that to this order the inequality is simply 

v 2 d d - 1 x > , (47) 



and is always satisfied. This result can also be seen by noting that the viscous order correc- 
tion to the energy ( 133|) is of order e 3 . Alternatively, we might have anticipated the above 
result on general ground. At order e 2 the only scalar quantities that we can build are v 2 , P 
and diVi. The pressure P appears in the NS equations via a derivative, thus, we can expect 
that it will not appear by itself, while d^Vi = 0. 

A non-violation of the inequality is only a necessary condition for not having finite-time 
blowups, still is it curious that it holds so simply for all incompressible flows. On the gravity 
side, we may expect that the gravitational backgrounds dual to the incompressible Navier- 
Stokes equations will not develop a naked singularity. Note also in comparison, that for 
the relativistic flows which are compressible the inequality is not obviously satisfied. This 
presumably suggests that singularities are more likely to be expected in compressible fluid 
motions. 

Finite-time blowups of three-dimensional nonrelativistic incompressible flows have been 
demonstrated with pure imaginary vector fields in [22] . It is amusing that in these cases the 
Penrose inequality is violated, unless we consider the absolute value of the flow field. 

C. Second order viscous flows 

Let us consider the first subleading correction in the 1/c expansion to the inequality of 
the incompressible case, which is now at order e 4 . The fluid flows are now compressible. For 
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simplicity we will analyze the black brane hydrodynamics and use ( |4"2]) . The energy reads 
„ d-1 f ( dv 2 ,„ dv 4 d 2 Pv 2 (d + l)b 2n A . , 



16vr^ J V d-1 d-1 d-1 (d-1) 2 , 

Note that at this order we cannot use the incompressibility condition and the last term does 
not vanish. The new transport coefficient r w does not appear at this order. 
The entropy in the zeroth order of derivatives is 

The entropy is corrected at second order in derivatives (1551) 

/ (A l( r tj (r lj + A 2 u ij u ij + A 3 R + B x D i( j Qi + B 2 DiU 0i ) d d ~ l x . (50) 

J 46 

The inequality reads 

5v 4 dPv 2 b (d+l) 2o 4 , x ^ 
— + — — + ^— — J -v 2 diV l d d x x > 
8 2 a 

hi J (^±^)d iVj dV + {^ L Y^-)d i v j d j v l - ^Yf^fd^x , (51) 
and is not obviously satisfied. 



V. DISCUSSION 

In this paper we derived an inequality, analogous to the Penrose inequality, for relativistic 
and nonrelativistic fluid flows in a general dimension. The inequality provides means to study 
singular flows. As in gravity, a non-violation of the inequality is only a necessary condition 
for the absence of singularities. Indeed, we saw that the inequality is always satisfied at the 
ideal fluid order, and in particular is not violated by relativistic shock waves. 

At the leading viscous order, the inequality may be violated by relativistic fluid flows, A 
construction of an example that violates the inequality or a proof that it is always satisfied 
is still an open problem. We saw that the inequality at this order is always satisfied by 
nonrelativistic incompressible flows. Hence, on the gravity side, we may expect that the 
gravitational backgrounds dual to the incompressible Navier-Stokes equations will not de- 
velop a naked singularity. The inequality may be violated at the next to leading viscous 
order by both relativistic and nonrelativistic compressible flows. This presumably suggests 
that singularities are more likely to be expected in compressible fluid motions. 
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The derivation of the inequality is rather general and can be applied to any closed dy- 
namical system that relaxes to equilibrium at late times. In the context of fluid dynamics it 
can be easily applied also to charged hydrodynamics. An interesting arena for a future study 
of the inequality is time dependent nonequilibrium relativistic dynamics. In this case the 
inequality may depend crucially on the definition of the entropy, which is not uniquely de- 
fined away from equilibrium (see e.g. 23] ). Finally, there are other interesting cases where 
one can consider the Penrose inequality, e.g. black hole horizons with more complicated 
topology. 
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